arXiv:astro-ph/0211505vl 22 Nov 2002 


Sensitive Limits on the Water Abundance in Cold Low Mass 

Molecular Cores 

Edwin A. Bergin 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

eberginOcfa.harvard.edu 
Ronald L. Snell 

Department of Astronomy, University of Massachusetts, Amherst, MA 01003 

snellQastro.umass.edu 

ABSTRACT 

We present SWAS observations of water vapor in two cold star-less clouds, 

B68 and Core D in p Ophiuchus. Sensitive non-detections of the lio — loi transi¬ 
tion of 0 -H 2 O are reported for each source. Both molecular cores have been pre¬ 
viously examined by detailed observations that have characterized the physical 
structure. Using these rather well defined physical properties and a Monte-Carlo 
radiation transfer model we have removed one of the largest uncertainties from 
the abundance calculation and set the lowest water abundance limit to date in 
cold low-mass molecular cores. These limits are x(o-H20) < 3 x 10“® (relative 
to H 2 ) and x(o-H20) < 8 x 10“® in B68 and p Oph D, respectively. Such low 
abundances confirm the general lack of ortho-water vapor in cold (T < 20 K) 
cores. Provided that the ortho/para ratio of water is not near zero, these limits 
are well below theoretical predictions and appear to support the suggestion that 
most of the water in dense low-mass cores is frozen onto the surfaces of cold dust 
grains. 

Subject headings: astrochemistry - ISM: molecules - stars: formation - submil¬ 
limeter 


1. Introduction 

Water is an important molecule in the interstellar medium (ISM) because it links chem¬ 
istry in the ISM to comets and planetary systems and provides crucial aid in constraining 


the chemistry of astrophysical systems (Bergin et ah 2000). Recently NASA’s Submillimeter 
Wave Astronomy Satellite (SWAS) and ESA’s Infrared Space Observatory (ISO) have pro¬ 
vided the hrst unambiguous glimpse of water in the ISM. Because the SWAS instrument is 
pre-tuned to the ground state rotational transition of ortho-H 20 , it is sensitive primarily to 
water in the quiescent or coldest regions of molecular clouds. In contrast, ISO sampled higher 
energy transitions that consequently probed warm gas found in shocks or in close proximity 
to young embedded stars. One of the surprises from the analysis of SWAS results is that the 
derived water abundance is unexpectedly low, with values x(o-Il 20 ) ~ 1 —10 x 10“® (relative 
to H 2 ) found in a variety of regions (Snell et ah 2000a,b,c). These abundances conflicted 
with theoretical expectations and several suggestions have been provided to account for the 
discrepancy, with the principal solution for cold gas being the freeze-out of water onto grain 
surfaces (Bergin et ah 2000; Spaans & van Dishoeck 2001; Viti et al. 2001; Charnley, Rodgers, 
& Ehrenfreund 2001). Moreover, ISO and SWAS detections of water vapor in absorption 
produced results in disagreement with the emission line analysis. Absorption lines are more 
straightforward to analyze and provide water abundances of a:(o-Il 20 ) ~ 10“®, in agreement 
with theory (Cernicharo et al. 1997; Neufeld et al. 2000; Moneti, Cernicharo, & Pardo 2001). 
These differences can be reconciled provided that the freeze-out of water is greater in denser 
regions, which are seen in emission, as opposed to the low density gas, with correspondingly 
longer depletion timescales, seen in absorption (Bergin et al. 2000). 

Nonetheless there are some well discussed uncertainties in the analysis of water emission 
lines detected by SWAS (Snell et al. 2000a,b,c,). These arise primarily from the fact that 
only a small column of water is required to produce optically thick emission. The situation 
is slightly simplihed, because the critical density for the lio — loi transition (~ 10® cm“®) 
is higher than typical densities found in molecular cores. Under this condition every photon 
that is collisionally created eventually escapes and the water emission, although optically 
thick, is effectively thin (Wannier et al. 1991; Snell et al. 2000a). This fact aids in constrain¬ 
ing the abundance, but crucial assumptions must be made regarding the source physical 
structure. Assuming that the lio — loi emission arises predominantly from high density 
gas, simple single component models were constructed (i.e. single density and temperature 
along the line of sight), resulting in the above abundance estimates (Snell et al. 2000a,b,c). 
However, uncertainties in the density are directly reflected in the abundance, and a single 
density characterizing the entire line of sight is itself an approximation. 

More realistic physical structures were examined for a few molecular cores using a Monte- 
Carlo radiative transfer code (Ashby et al. 2000a). This analysis supported the simple single 
component results. However, there are complications. For instance in S140, a separate anal¬ 
ysis argues that a clumpy cloud model is more appropriate, although the resulting abundance 
estimates are similar (Spaans & van Dishoeck 2001). Sources that offer the most promise to 
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set good limits or determinations of the water abundance are therefore ones that are largely 
isolated, and have well established physical properties. In general, giant molecular cloud 
cores, such as S140, do not £t this criteria as they have been suggested to have clumpy 
physical structures which are difficult to model definitively. On the other hand, advances 
in the resolution and sensitivity of infrared (IR) and submillimeter continuum arrays have 
begun to provide a wealth of information on the physical structure of low mass cores. Studies 
using these instruments effectively apply similar techniques in which azimuthally averaged 
dust column density profiles are combined with assumptions regarding the geometry and a 
gas-to-dust ratio to determine the H 2 density profile. 

We present here SWAS observations of two molecular cores, B 68 and p Oph D, that 
have been subject to these techniques and thus have rather reliable estimates of the source 
physical structure. Using these physical models and SWAS observations we have set the 
lowest limit to date on the water abundance in cold low mass cores. The limits are below 
pure gas-phase chemical predictions and support the assertion that most of the water is likely 
frozen onto grain surfaces in cold dense cores of molecular clouds. 


2. Observations and Resnlts 

Between Feb. 25, 2001 to Mar. 8 , 2001 SWAS observed B 68 in the lio — loi transition 
of 0 -H 2 O at 556.936 GHz for a total of 28 hours of on-source integration. Similarly during 
Sept. 7-10, 2001 SWAS observed p Oph D for a total of 24.6 hrs (on-source). System 
temperatures were typically 2500 K with minimal scatter around that value. The spacecraft 
was used in nod mode, involving alternately nodding the spacecraft to an off-source position 
free of emission. For B 68 the elliptical 3l3 x 415 beam was centered on a = 17:22:38.2, 5 = 
-23:49:34 (J2000) and towards p Oph D, a = 16:28:30.4, 5 = -24:18:29 (J2000). All data 
were reduced with the standard SWAS pipeline described by Melnick et al. (2000). Towards 
both cores, along with 0 -H 2 O, SWAS simultaneously observed transitions of [C I] (^Pi —> 
^Pq), ^^CO (J = 5 ^ 4), and O 2 (83 —I 2 ). In this work we primarily use the 0 -H 2 O and 
^^CO data, with a velocity resolution of 1.0 kms“^, and a sampling of 0.6 kms“^. All data 
are presented here on the T^ scale and for subsequent analysis are scaled by the main beam 
efficiency of 90% (Melnick et al. 2000). 

Figure 1 shows spectra of the lio — loi transition of 0 -H 2 O and J = 5 — 4 taken 
towards B 68 and p Oph D. It is apparent that there are no water detections toward either 
source. In B 68 the 3cr upper limit is T^ = 36 mK and in p Oph D the limit is 45 mK. The 
J = 5 — 4 transition of ^^CO was not detected towards B 68 with a similar limit as 0 -H 2 O, 
confirming the cold (T < 10) nature of this source. ^^CO J = 5 — 4 is detected towards p 
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Oph D with an integrated intensity of 0.30 K km s a peak antenna temperatnre of = 
150 mK, and Av = 1.86 km s“^. 


3. Monte-Carlo Models of Water Emission 

To estimate the water abnndance we draw npon previous determinations of the density 
structure derived from observations of dust in emission or absorption (Motte, Andre, & Neri 
1998; Alves, Lada, & Lada 2001; Bacmann et ah 2000). Towards both clouds, the SWAS 
beam encompasses the entire area seen in transitions of molecules such as C^®0, CS, and 
N 2 H+. Higher spatial resolution observations of molecular cores in these other tracers are 
able to determine the extent of emission, which along with the radial density prohle with 
radius, permits a more accurate analysis of abundances (Jessop & Ward-Thompson 2001; 
Tafalla et ah 2002; Bergin et al. 2002; Hotzel et ah 2002). Because of the large beamsize our 
initial procedure is to assume a simple model with an abundance of 0 -H 2 O that is constant 
with cloud depth. However, chemical models of centrally condensed objects predict that 
lower abundances should exist in the dense central regions where molecules freeze onto grain 
surfaces more frequently (Rawlings, Hartquist, Menten, & Williams 1992). To account for 
this possibility we have also examined whether a water abundance prohle predicted by theory 
is in conhict with our observations. This model is applied to B68 which has a large amount 
of additional observational constraints. 

The water abundance prohle, along with the adopted density, temperature, and veloc¬ 
ity width prohles, is used as input into the spherical one-dimensional Monte-Carlo radiative 
transfer code discussed by Ashby et al. (2000b).^ The adopted physical prohles for each 
source are motivated below, but in each case we have assumed a static cloud. The radia¬ 
tive transfer model determines the expected emission spectra, which is compared to the 
observed 3a upper limit on the antenna temperature. The water abundance is then iterated 
until the predicted emission matches the limit. To account for oversampling in the SWAS 
spectrometer, we have convolved the model data with a Gaussian that has a width of 1.45 
MHz. Because of the weak IR continuum radiation held in both cores the ehects of infrared 
pumping are negligible. 


^ Rates of excitation of water with ortho-H 2 are often an order of magnitude greater than rates for lower 
energy para-H 2 (Phillips, Maluendes, & Green 1996). For our calculations of cold clouds we have assumed 
an ortho/para-H 2 ratio of 0.1 and extrapolated the rates down to 10 K. If the ratio is lower our results will 
be essentially unchanged; a higher ratio would result in lower abundance limits. 
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3.1. B68 

The B68 globule is located at a distance of 125 pc (de Gens, de Zeeuw, & hub 1989). 
Alves, Lada, & Lada (2001) examined this isolated pre-stellar core using near-infrared extinc¬ 
tion techniques. The resulting extinction map was used to constrain the radial density prohle, 
which is well £t by a pressure confined, self-gravitating cloud near equilibrium (Bonnor-Ebert 
sphere). Recently, Bergin et al. (2002, hereafter BAHL02) and Hotzel et al. (2002) have ex¬ 
amined this cloud in the low-J mm-wave transitions of various molecules. Both found that 
in the dense core center, molecules are systematically depleted onto grain surfaces. 

We use the model described by BAHL02 to estimate the water abundance. In this model 
the column density is constrained by the visual extinction data and the line width by 
emission (BAHL02). The assumed density profile is that derived from near-IR extinction 
measurements, while the temperature structure is that for dust in a Bonnor-Ebert sphere 
derived by Zucconi, Walmsley, & Galli (2001).^ The density and temperature structure are 
provided in Figure 2. This model assumes a static cloud, but includes contributions from 
thermal and turbulent line widths with the latter increasing as a function of radius. 

As a check on the assumptions regarding the cloud physical structure we first examine 
the non-detection of ^^GO J = 5 — 4. To model this emission we use the depleted G^®0 
abundance profile given by BAHL02 scaled by a factor of 7.8 (accounting for ^^G/^^G = 75; 
leQ/iSQ = 500). Placing this into the radiative transfer code produces a peak temperature of 
T^ ~ 40 mK, which is close to the observed 3cr limit of 36 mK. Since the ^^GO analysis does 
not significantly conflict with current limits, we then use the iterative procedure outlined 
previously to estimate the 3cr limit on the water abundance in B68, which is a;(o-H20) 
< 3 X 10“®. The resulting emission spectra from the radiative transfer model are shown in 
Figure 1.^ 

BAHL02 linked a radiative transfer code to a chemical model including the effects of 
molecular depletion to set limits on the freeze-out of G^®0 molecules. Water is more tightly 
bound to grain surfaces when compared to GO (Sandford & Allamandola 1988). Therefore 
it is more difficult to remove and can be expected to be significantly depleted in the center, 
perhaps with larger depletions than GO. To see if current chemical theory is in in conflict 
with our observational results we compare the predictions of a gas-grain chemical model to 


^However, for the temperature, a uniform reduction of 2 K was required to match the multi-transitional 
N 2 H+ observations of BAHL02. 

^If the core radius is allowed to extend beyond 0.06 pc, as estimated by the near-IR analysis, then the 
abundance limit will be reduced. 
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the observations. 

To this end we use the gas-grain chemical model discussed in BAHL02. The only 
difference in the chemistry is the additional inclusion of the grain surface formation of water 
via hydrogenation of oxygen as required by previous studies (Bergin et ah 2000; Viti et al. 
2001; Charnley, Rodgers, & Ehrenfreund 2001). This model incorporates the cloud density 
and temperature prohle shown in Figure 2. We also include the observed cloud extinction 
prohle which increases to Ay = 17™ at the core center. The chemical model predicts the 
prohle of water abundance with cloud radius which, along with the same physical prohle, 
are placed as inputs to the radiative transfer calculations. Due to the inclusion of gas-grain 
interactions the chemical model does not reach a steady state and the predicted abundance 
prohle strongly varies with time. Fortunately BAHL02 hnd that the observed and N 2 H+ 
emission in B68 can be simultaneously reproduced by the model at t ~ 7.6 X 10"^ yr (later 
times predict increasing depletion and, in consequence, less emission than observed). 

We therefore adopt this time in our analysis to examine whether the H 2 O abundance prohle 
predicted by the same model is in conhict with SWAS observations. 

In Figure 3 we provide the predicted water vapor and ice abundances as a function of 
depth. Also shown are the abundance of atomic oxygen and the best ht constant abundance. 
The chemical model predicts that the water vapor abundance sharply declines at the cloud 
edge due to photodissociation, rises to a peak at Ay ~ 2™, and then uniformly declines 
towards higher depths. As seen in Figure 3, nearly all of the water and oxygen (in the 
form of H 2 O) is frozen onto dust grains. Placing the results of this model into the radiation 
transfer code we predict a peak temperature in the 0 -H 2 O 557 GHz line of = 18 mK 
below the 3a upper limit. As long as the age is > 10^ yr our observations are consistent 
with this model. Thus current predictions of water abundance prohles by chemical models 
are not in conflict with observations. Given the potential of B68 as a template for testing 
chemical models it will be useful in the future to re-examine these results in light of additional 
molecular observations. 


3.2. p Oph D 

The p Oph D molecular core is a star-less object located in the Ophiuchus complex 
at a distance of 160 pc. The core was detected in absorption against the galactic mid- 
IR background using ISOGAM (Bacmann et al. 2000) and in 1.3 mm continuum emission 
(Motte, Andre, & Neri 1998). These studies have provided good constraints on the core 
density prohle, although there are differences in the details. Both agree that the density 
declines as p{r) oc beyond r ~ 3300 — 4000 AU; within this radius the density is found 
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to be constant, ranging from ricst = 3 x 10^ cm“^ (Bacmann et al. 2000) to ricst = 9 x 10^ 
cm“^ (Motte et al. 1998). Both studies derive a core radius of ~ 13,000 AU (0.063 pc). 
However Motte, Andre, & Neri (1998) find a sharp edge to the South-West but, in the East- 
West direction, the dust emission merges with the ambient cloud. With these differences core 
mass estimates range from 2-5 M©. In the following we examine both density distributions. 

Unlike B68, where a detailed model exists, there is no information regarding the radial 
dependence of the velocity held or temperature in p Oph D. For temperature, the sub¬ 
millimeter dust continuum survey of Ophiuchus by Johnstone et al. (2000) did not include 
core D, but for numerous other cores in the cloud they hnd Tdust ~ 10 — 30 K, which may 
be considered the expected range. Given the lack of detailed information, we adopt the 
expression for the equilibrium dust temperature given by Burton, Hollenbach, & Tielens 
(1990), further assuming that "Tdust = Tgas- This expression provides an estimate of the dust 
temperature depending on the extinction and of the local enhancement of the ultra-violet 
radiation held, Gq. For the Ay(r) prohle we integrate the observed density prohle with 
radius in a pencil beam along the line of sight. Because of the presence of early type stars 
the local radiation held is enhanced throughout the p Oph cloud and detailed modeling of 
the far-infrared emission by Liseau et al. (1999) hnds Gq = 20 towards core D. 

The procedure is to hrst examine how well the physical model(s) reproduce the SWAS 
^^GO J = 5 — 4 emission and then apply the “best” model to the water emission. If we use 
the Bacmann et al. (2000) density prohle, the above temperature structure, and a constant 
turbulent width of 1.5 km s“^ we can match the observed ^^GO emission provided a;(^^GO) 
~ 3 X 10“®. This abundance is higher than the expected range for undepleted gas.^ Since 
this source is a cold pre-stellar object we expect to see some evidence for ^^GO depletion. For 
instance in B68, and several other similar pre-stellar objects, large depletions are observed 
(Bergin et al. 2002; Jessop & Ward-Thompson 2001; Tafalla et al. 2002). 

Adopting the denser prohle of Motte, Andre, & Neri (1998) we can reproduce observed 
^^GO emission at a lower abundance of ~ 4 xl0“^. The predicted emission spectra from 
this model is shown superposed on the observational data in Figure 1. For reproducing the 
water emission, the Motte, Andre, & Neri (1998) model is favored due to the anticipation 
of GO depletion. Adopting this model for the lio — loi observations we derive a;(o-H20) 
< 8 X 10“®, relative to H2(3cr). The model that matches the 3a limit is also provided in 
Figure 1. If we use the Bacmann model the abundance limit is a:(o-H20) < 1 x 10“^. 


"‘Using the C‘®0 abundance derived in the extended Ophiuchus cloud by Frerking, Langer, & Wilson 
(1982) with plausible isotope ratios (‘®0/‘®0 = 500; = 45 - 90) the abundance is expected 

to be 1 - 2 X 10“®. 



4. Conclusions 


We have derived upper limits to the water abundance in two clouds with well described 
physical properties. These abundance limits are, a;(o-H20) < 8 x 10“® for core D in p Oph 
and x(o-H20) < 3 x 10“® in B68. These limits are below those previously set for a cold star¬ 
less object by Snell et ah (2000c) (TMC-1; a;(o-H20) < 7 x 10“®). In this Letter two cores 
have been subject to a more detailed and careful study which essentially conhrms the general 
lack of ortho-water vapor in cold (T < 20 K) molecular cores. Provided the ortho/para ratio 
of water is > 0.03 then these results are in agreement with the assertion that the water 
abundance in low mass objects is well below the predictions of pure gas phase chemistry. In 
the case of B68, these observations are compared to theoretical predictions of a gas-grain 
chemical model, which have been directly placed into the radiative transfer calculations. We 
hnd that these model predictions are not in conflict with observations. However, due to the 
large SWAS beam, which encompasses the entire extent of the observed molecular emission 
seen in other tracers, we cannot discriminate between simple constant abundance models 
and those with complex abundance structure constrained by theory. Given the wide-spread 
molecular depletion found in B68 and similar cores (Bergin et ah 2002; Hotzel et ah 2002; 
Tafalla et ah 2002), it is likely that most of the water in low mass dense cores is frozen on 
the surfaces of cold dust grains. 


We acknowledge Gary Melnick for a thorough reading of the manuscript. E.A.B. is 
grateful for the help and collaboration with Gharlie Lada, Joao Alves, and Tracy Huard 
which greatly aided the B68 analysis. For these data we are grateful to entire SWAS team 
and acknowledge support from NASA’s SWAS Grant NAS5-30702. 
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Fig. 1.— Spectra of the lio — loi transition of 0 -H 2 O and J = 5 — 4 towards the B68 
dark clond (left) and p Oph core D (right) shown as a thick solid line. The thin solid line is 
the resnlt of the excitation model described in §3.1 and §3.2, while the dashed line denotes 
the systemic velocity. The model is shifted to the proper source velocity. J = 5 — 4 

is detected in p Oph D, but all other observations are non-detections. In the hgure a linear 
baseline has been subtracted from the observational data. The model spectra are shown 
with continuum included which, due to the small level, has not been subtracted. 
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Fig. 2.— Adopted density and temperature profile with radius in B68 (left) and p Oph D 
(right) molecular cores. In the latter core we present the density prohle derived by Motte et 
ah (1998) and Bacmann et al. (2000). 



Fig. 3.— Theoretical predictions of molecular abundance as a function of cloud depth (ex¬ 
tinction) from a gas-grain chemical model. Also shown is the best £t constant abundance 
model for the B68 cloud. Details of this model are given in §3.1. 
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